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' We theoretically investigate electron transport through an Aharonov-Bohm interferometer con- 

taining laterally coupled double quantum dots. We introduce the indirect coupling parameter a, 
which characterizes the strength of the coupling via the reservoirs between two quantum dots. 
|a| = 1 indicates the strongest coupling, where only a single mode contributes to the transport in 
the system. Two conduction modes exist in a system where \a\ ^ 1. The interference effects such 
CO , as the Fano resonance and the Aharonov-Bohm oscillation are suppressed as the absolute value of 

the parameter a decreases from 1. The linear conductance does not depend on the flux when a = 
since it corresponds to independent coupling of the dots to the reservoir modes. 
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PACS numbers: 73.23.-b, 73.63.Kv, 73.40.Gk 



I. INTRODUCTION 



Quantum phase coherence in mesoscopic systems has attracted the attention of many physicists. Quantum phase 
coherence is detectable by interference experiments. In interference experiments with an Aharonov-Bohm (AB) ring 
containing a quantum dot (QD), quasi-periodic modulation of the tunneling current as a function of the magnetic 
' flux through the ring has been experimentally demonstratediiSi^. This reflects the fact that the phase coherence is 
P5 ■ maintained during the tunneling process through a QD. The Fano effect^ is another important interference effect in 
O mesoscopic physics. The Fano effect occurs in a system in which discrete and continuum energy states coexist^. 
, ^ , , Recently the AB oscillations of a tunneling current passing through a laterally coupled double quantum dot (DQD) 
system were observed by Holleitner et al. using a lateral DQD- . Such oscillations were also observed by a group 
^sj , including two of the present authors using a vertical DQDi in a weak inter-dot tunnel coupling regime. Motivated 
J> ' by these experimental results, electron transport through such a system has been investigated theoretically^iSiiS. 
Moreover, such systems are very interesting since it has been theoretically predicted that cotunneling currents passing 
■ through spin-singlet and triplet states have different AB oscillation phasesii. DQD has been attracting attention as 
[ an important device structure for entangled spin qubit operationsiSiiii^. 

. In this paper, we consider the transport through an AB interferometer containing a laterally coupled DQD. Although 
\Q ' the electron spin and interaction effects are crucial in the previous theoretical proposals, here we disregard them and 
, focus on the single particle interference properties. Instead, we introduce the indirect coupling parameter a, which 
characterizes the strength of the coupling via the reservoirs between two QDsi^. A system with a maximum coupling 
^ ] |a| = 1 has already been widely studied theoreticall}iSi2ii£. In actual systems, however, such a case is very special 
and most experimental situations correspond to |q;| < 1. We found that the number of conduction modes is two 



' I except when \a\ = 1, where only a single mode contributes to the transport. The situation where a = has also 



been explored in the context of the orbital Kondo problemiSiii. We calculate the tunneling current through the 
Q DQD systems in terms of Green's function techniques for non-interacting systemsiSiiS. By using an exact expression 
^ , of the tunneling current, we examine AB oscillations through DQD from the weak inter-dot tunnel coupling regime 
to the strong inter-dot tunnel coupling regime. The visibility of AB oscillation (i.e., the ratio of the AB oscillation 
amplitude and the current maximum) decreases as the absolute value of the indirect coupling parameter a decreases. 
The current is independent of the flux at a = 0, and the visibility becomes zero. As a function of energy, the linear 



J-j [ conductance indicates an asymmetric Fano lineshape when \a\ = 1, when the energy of the (anti-)symmetric state 
is near the Fermi level of the reservoirs. However, the Fano resonance is suppressed as |a| decreases and finally the 
linear conductance exhibits two peaks at a = corresponding to the Breit-Wigner resonances by the symmetric and 
anti-symmetric states. 

This paper is organized as follows. In Sec. m a standard tunneling Hamiltonian is employed to describe an AB 
interferometer containing a laterally coupled DQD. We introduce the indirect coupling parameter a in Sec. IIII Al We 
derive the linear conductance at zero temperature in Sec. IIII Bl In Sec. IIII CI we provide the energy dependences of 
the linear conductance at zero temperature and we show that the Fano resonance is suppressed in a system where 
\a\ 1. In Sec. IIII Dl we discuss the AB oscillations of the linear conductance at zero temperature. In Sec. IIVI we 
extend the argument where the magnetic flux and tunnel coupling between reservoirs and QDs are asymmetric. All our 
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results are summarized in Sec. In Appendix^ we give the detailed derivation of the indirect coupling parameter 
a for the two-dimensional and three-dimensional reservoirs. In Appendix^ by diagonalizing the transmission matrix, 
we show that there are two conduction modes when \a\ ^ 1. However, there is only a single conduction mode when 
|a| = 1. When a = 0, the two conduction modes correspond to the Breit-Wigner resonances by the symmetric and 
anti-symmetric states. 



II. MODEL 



We consider an AB interferometer containing a DQD coupled to two reservoirs as shown in Fig. Only a single 
energy level in each QD is assumed to be relevant and we ignore the spin degree of freedom. We model this system 
with the following Hamiltonian: 

H = Hji + Hdqd + Ht, (1) 
where Hji describes Fermi seas of non-interacting electrons in two reservoirs 

Hr= ^ '^'ivkCuk^ Cuk- (2) 

ve{ux} k 

Here t^^. is the energy of conduction electrons with wave number k in the reservoirs v = U, L and the operator c^k 
(cyfc^) annihilates (creates) an electron in reservoirs. Hjjqd is the Hamiltonian of the isolated DQD 
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Hdqd = ^ ^odi^di -I- tc {di^d2 + h.c^ . (3) 

i=l 

Here eg is the single-particle energy level of two QDs2£ and di {di^) annihilates (creates) an electron in the ith QD 
{i = 1,2). The second term represents direct tunneling between two dots. The energy levels of the dots split to form 
symmetric and anti-symmetric states because of inter-dot tunnel coupling, t,.. Their energy levels are Cs = cq + tc and 
ta = ^0 — tc, where we choose the gauge such that is real and negative. Ht describes the tunneling Hamiltonian 
between reservoirs and QDs 
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\t^^le^'^cuk'd,+t^^le-^'^cukU, 



+trk^'^^Lk^d2 + t^Lle-'^CLk^di + h. 

2 



iye{U.L} k 1=1 

where t^^l is real. The factors e^*^^ indicate the effect of the magnetic flux {(j)^ = 27r<I'y/$o is an AB phase in each 
subcircuit, where is the magnetic flux threading through each subcircuit as shown in Fig. 1, and = h/e is the 
magnetic flux quantum). 



III. TUNNELING CURRENT 



For non-interacting systems, the tunneling current can be written asi^ 

I = tJ delfuie)-fLie)]Tr{G^ie)T^ie)G^ie)r^ie)} 

= Ij de[fu{e)~fL{e)]T{e). (5) 

Here /i/(e) = 1/ [l + e'^^^"-'/'^^"^] with v = U/L is the Fermi-Dirac distribution function, being the chemical 
potential of the reservoir v. Boldface notations indicate 2x2 matrices, where G'^^°'\t) is the retarded (advanced) 
Green's function of a DQD, r'^^^' is the linewidth function defined by 

rr,.(6) ^ 2^^t[;f(0,)t[/2(0.)'5(e - e^k), (6) 

k 

and T(e) is the transmission probability at energy e. 
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FIG. 1; Schematic diagram of an AB interferometer containing a laterally coupled DQD. The magnetic fluxes threading the 
upper and lower subcircuits are <^u and respectively, and cause the AB effect, s is the distance between the two QDs. 
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FIG. 2: Model to estimate the tunneling coupling between the reservoirs and the two dots. S indicates the area in the 
tunneling potential. The radius of the QD is Ro and the center position of ith dot is {xi,yi,0). 



A. Indirect coupling parameter a 

Here wc estimate the linewidth function based on Ref. |^ and focus on the three-dimensional (3D) reservoirs. The 
detailed derivation of the indirect coupling parameter a for not only 3D but also two-dimensional (2D) reservoirs is 

given in Appendix ^ According to the Bardeen's formula™, the flux independent tunneling matrix element t^jl is 
given by 

where ^ is the effective mass, the 2D area 5* for the integration is in the tunneling potential, (puki^) is the wave 
function of evanescent mode of the reservoir i/, and <Pd\''^) is the wave function of a localized electron in the dot i. 

We calculate the tunneling matrix elements for the model shown in Fig. |21 where we consider only the upper 
reservoir. As shown in Fig. |3 the boundary of the reservoir is a flat surface, the dot size is Ro, and the potential 
profile is 

{0 (in the reservoir, y < 0) 

EtiUdotir-n) i\f^n\<Ro) , (8) 

U (others including the area 5") 

where fi = (x^, 0) is the center of the ith dot. Then the wave functions of reservoir modes are 

cl^Ukir) - -^6^'=^-+^'=^^ sin e--^ (9) 
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where the decay constant k, is 



and the scattering phase shift 9 satisfies sinfl = tJ+k ^ • ^ system size and fc||^ = /cx^ + kz^. For the dot, we 

consider only the symmetric mode, namely zero angular momentum state, 

0|,*)(r>i?o) = Z?3^^e--l^~-^~"l (11) 
a\r-ri\ 

where + ^||^ = '^^ ^^'^ is a constant related to the detail of the potential inside the dot Udot- Using Eq. lO, 
the tunneling matrix elements are 



i'ui = :T77W«in - ^e-^-+--. (12) 



Then we calculate the linewidth function r^(e) using Eq. The diagonal elements r^(e) are real, which are 

equal since we assumed yi = y2 = TJd- We define the indirect coupling parameter as follows: 

1-2 I, 2 

We found the indirect coupling parameter at the Fermi energy /i^/ = with t/ ^ fiu^ 

a^itiu)- ^Jiikps), (14) 
where ji is the first order spherical Bessel function and for U > fijj 

" - (,2 + (2,^)2)3/2' (15) 

where s = |xi — a;2|. We find \oP\ < 1 and decreases with s. The oscillation with kps shown in Eq. (|14|l originates 
from the assumptions of ballistic transport in the reservoirs and the flat reservoir surface. This oscillation becomes 
weaker and a decreases monotonically with s if we consider curved or rough boundaries. We expect that the above 
discussions can be extended into the general model as shown in Fig. ^ We restrict our discussion to a symmetrically 
coupled DQD system, r^j(e) = r22(e) = 7iy(e) for simplicity. Then, using a basis of localized states, we use the 
linewidth functions 

in the following argument. 



B. Linear Conductance 



We obtain the following retarded Green's function using the equation of motion, 

1 -1 



(e-eo)/n + i(7c/(e)+7L(e))/2 + i(a^(e)7c/(e)e-*^^ + a^(e)7L(e)e'^^ )/2 

-t,/h + i(a^(e)7c/(e)e'^^' + a^(e)7L(e)e-''^^)/2 (e - eo)/h + i{ju{(-) + 7L(e))/2 



(17) 



where is the retarded Green's function of an isolated DQD. The advanced Green's function is G°'{e) = [^"(e)]^. 
The linear conductance at zero temperature is given by the Green's functions and the linewidth functions at the Fermi 
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energy (fiu = fi^). In the following, we choose the Fermi energy as the origin of energy. To clarify the physical image, 



= a and 



we consider the symmetric situation where 7c/ (0) = 7_l(0) = 7, a 
IIVI Then, we have the following transmission probability at the Fermi energy: 

T(0) = Tr{G''(0)r^(0)G''(0)r^(0)} 

2(^17)2 eo^{l + o? COS0) - 4aeotc cos + (1 + or-)tc^ + (h-^fil - o?) | 



/2 except for in Sec. 



1 — cos"' 



e 2 +2 



(^7)2 { 



1 — cos2 



I)}] +4(^7)^ 



eo — aic cos ( f 



This is the main result of this paper. Then, the linear conductance is given by 

G(eo,ic,«» = Ggr(0), 



(18) 



(19) 



where Gq = 



^'jh. This linear conductance has the symmetry of G(— eo,tc,Q!, 



G(eo, tc, a, + 27r). Therefore, in 



the following, we consider only eo > since we can obtain the results for ep < with adding 2ti to the AB phase of 
the linear conductance for eo > 0. From the flux dependence of the conductance, we find that the phases of the AB 
oscillations change by I-k via the resonance except when a = 0. Similarly, we consider only a > since the linear 
conductance has the symmetry of G(eo, — o;, (/)) = G(ep, ic, o;, + 27r). In particular, Eq. H19|) provides the same 
result obtained in Ref. [^j when a = 1 and that in Ref. llSl when = 0. 



C. Energy dependence 

We define the total magnetic flux threading through the AB interferometer as $ = In Fig. |31 we show the 

energy eo and the inter-dot tunnel coupling tc dependence of the linear conductance for a fixed flux $/<I>o = 1/8. 
While we use this flux value as an example, the following discussions are qualitatively valid except when $/$o is an 
integer which is discussed separately. Here we define the weak coupling regime with the condition |eo| > \tc\ and the 
strong coupling regime with the condition |eo| < \tc\- The former corresponds to regions I and III in Fig. El^c), and 
the latter corresponds to region II in Fig. Etc). When a = 1 (see Fig. E^a)), there is only one mode that contributes 
to the conduction (see Appendix IbI) and the conditions exist for a perfect transmission G/Gq = 1 whose curve has 
the shape of hyperbolas eo^ — tc^ = —{fvy)"^ sin^(0/2) (black broken curve in Fig. E^a)) in the strong coupling regime 
and of perfect reflection G/Gq = whose curve has the shape of cq cos((/)/2) — = (white broken line in Fig. El^a)) 
in the weak coupling regime. 

When a 7^ 1, effectively two modes contribute to the conduction (see Appendix 0. Thus, the maximum value of 
G/Gq may exceed 1 (see Fig. Elb), (c) and the changes in the scale). Moreover, at a = 0, the linear conductance 
becomes symmetric with respect to eo = as shown in Fig. EI^c) since two conduction modes have transport properties 
through the symmetric state (eo + 1^) and anti-symmetric state (eo — tc): respectively (see Eqs. (IB9|I and (|B10ll with 
setting e = 0). 

We inspect the eo-dependences of the linear conductance in detail. Figure 01 shows the eo-dependences of the linear 
conductance when tc/^7 = —1 and ^/^o = 1/8. When a = 1 (solid curve), the linear conductance has an asymmetric 
Fano lineshape near the resonance via the anti-symmetric state. The behaviors of the linear conductance at a = 1 
can be explained by the Breit-Wigner resonance via the symmetric state (eo = — tc) and the Fano resonance via the 
anti-symmetric state (eo = tc) as discussed in Ref. 0. The Fano effect is suppressed as a decreases and the independent 
Breit-Wigner resonances via the symmetric and anti-symmetric states realized at a = lead to two peak structures as 
shown in the dashed curve in Fig. 0] Although we considered tc/h-^ = —1, tc merely decides the resonance positions 
of the symmetric and anti-symmetric states and the above discussion is valid for any finite tc. In the following, we 
discuss the suppression of the Fano effect. When the energy of the anti-symmetric state is near the Fermi energy, the 
transmission probability at the Fermi energy shows the following generalized Fano forni24 

T(o)c.r,i!^i±^S±ii^^. (20) 

1 

Here the renormalized energy of the anti-symmetric state is defined by Ca = -^p-. We define the linewidth functions 
of the symmetric and anti-symmetric states, respectively, by 



r, =7<il + acos( IJ ^ r„=7jl-acos( I ) }, (21) 




~3 -2.5 O l-S -I -U.5 

(c) tjny 



FIG. 3: Gray-scale plots of the normalized linear conductance G/Gq for $/"l"o = 1/8. (a)a = 1. (b)a = 0.5. (c)a = 0. The 
gray dotted lines indicate the boundaries between the weak coupling regime (regions I and III, |eo| > \tc\) and strong coupling 
regime (region II, |eo| < |ic|). 



and the transmission probability of the background channel Tf, is the Breit-Wigner form via the symmetric state 
peaked at = + 2tc ^ '^tc alone 



(22) 
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FIG. 4: The eo-dependences of the hnear conductance when tc/h"f = —1 and "I>/<I>o = 1/8. The soUd, dotted, and broken 
curves indicate a — 1, 0.5, and 0, respectively. 
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FIG. 5: The suppression of Fano anti-resonance due to the indirect coupling parameter ct when tc/h-y = — 1 and ^/^o — 1/8. 
The solid, dotted, and broken curves indicate a — 1, 0.99, and 0.98, respectively. 



The Fano parameter q is positive and given by 



2a2sin2(|) 




1 — q;2 cos^ 




hVs 



(23) 



and the additional factor contributes to the transport only when a ^ 1 and is given by 

1 -\- a I 1 — 'T^ CDS li] I 



1 ~ COS" 



(I) 



2 ■ 2 

a sm 



(t) 



1 



(24) 



According to Eq. H20|l . when a = 1, the transmission probability at the Fermi energy provides the Fano lineshape 
and reproduces the result of Ref. 0. Moving away from a = 1, however, the transmission probability deviates from 
the Fano lineshape because of the additional term in Eq. (|2()() . Therefore, the Fano effect is suppressed and the Fano 
anti-resonance disappears as shown in detail in Fig. |S| For smaller \(j>\, the contribution of Ta is larger as seen in 
Eq. H24|l . Although here we discuss the situation of (f>/$o — 2ri, where n is an integer, the role of the symmetric and 
anti-symmetric states are interchanged when $/$o — (2n + 1) and the Fano parameter is negative and given by 



(25) 



2a^ sin^ 






\ — oP- cos^ 




fi-Fa 



Then we discuss the special situation that arises when $/<&o is an integer. First, we consider a = 1 where only a 
single conduction mode contributes to the transport. When $/$o = 2n electrons can only couple with the symmetric 
state since electrons tunnel from reservoirs into QDs with the same tunnel coupling under this condition as shown in 
Fig. EJa). Similarly, when $/$o = 2n + 1. electrons can only couple with the anti-symmetric state (sec Fig. Elb)). 
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(a) (b) 

FIG. 6: The tunnel coupling between reservoirs and QDs. The tunneling parameter t represents t^^l in Eq. (|7|l. (a) The tunnel 
coupling between reservoirs and QDs at "l>/$o = 2n. (b) The tunnel coupling between reservoirs and QDs at $/$o = 2n + 1. 
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FIG. 7: The linear conductance as a function of single electron energy level eo- The solid curve indicates the linear conductance 
for a = 1, tc/hj — ~1, and $/<I>o ~ 2n. The dotted, and broken curves indicate the linear conductances of two modes Ti and 
72 defined by Eqs. ijBl^ and (|B2j, respectively, for a = 0.8, tc/h"/ = —1, and $/$o = 2n. 

We consider the effects of a on the energy dependence of the linear conductance when <I>/<i>o = 2n. When a ~ 1, 
the linear conductance exhibits the lineshape of a Breit-Wigner type transmission via the symmetric state as shown 
in Fig. However, when a ^ I, the linear conductance has a side peak structure via the anti-symmetric state in 
addition to the Breit-Wigner like transmission via the symmetric state. This is because transport becomes possible 
by both the symmetric and the anti-symmetric states. The linewidth functions in the symmetric and anti-symmetric 
base are 

This is equivalent to Eq. I|21|) with (p ~ Ann. Thus, electrons can couple with both the symmetric and anti-symmetric 
states when a =/= 1. The ratio of the fuU-width-at-half-maximum of the two conductance peaks shown in Fig. [7|is the 
ratio of 1 -|- a and 1 — a. Similarly, when <&/$o = 2n + 1, the roles of the symmetric and anti-symmetric states are 
reversed. 



D. AB oscillations 

In this section, we discuss the magnetic flux dependence of the linear conductance. In particular, the linear 
conductance is independent of the flux at a = and the visibility of AB oscillation (i.e., the ratio of the AB 
oscillation amplitude and the current maximum) becomes zero in all situations. 

1. Zero coupling 

First, we consider the weak coupling limit, namely the absence of direct tunneling, tc = 0. AB oscillations are 
shown in Fig. |S1 When a = 1, the linear conductance shows a usual AB oscillation with a visibility of 1 and a period 
of 2tt because of the effect of the interference between the probability amplitude associated with one path through 
QD 1 and the probability amplitude associated with another path through QD 2. 
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FIG. 8: AB oscillations in the absence of inter-dot tunnel coupling. The solid, dotted, and broken curves indicate a — 1, 0.5, 
and 0, respectively, when eo/hj = 0.5. Inset: a dependence of the visibility of zero coupling. The solid, dotted, and broken 
curves indicate eo/h'y — 0.5, 1, and 2, respectively. 



/r 




FIG. 9: AB oscillations in the weak coupling regime, eo/h'y — 2 and tc/h'y — —1. The solid, dotted, and broken curves indicate 
a = 1, 0.5, and 0, respectively. 



For general a, the visibility v is given by 

(1 - a^) \{eo' - (^7)^(1 - a^)V + m)W 



[eo2 + (ft7)2] [eo2(l + a2) + (ft7)2(l 



(27) 



The indirect coupUng parameter dependences of the visibility is shown in inset to Fig. |S| The visibility decreases as 
a decreases. In particular, when |eo/?i7| 3> 1, — 2q;^/(1 + 0^). While, when |eo/?l7| <C 1, w — a^- For any cq, the 
visibility v is nionotonically increasing function of a. 



2. Weak coupling regime 

Here we consider the weak coupling regime, |eo| > \tc\. AB oscillations in this regime are shown in Fig. Elfor various 
parameters a. There are completely destructive interferences for a = 1. The visibility decreases as a decreases as 
found in the weak coupling limit. The AB oscillation period becomes 47r since the linear conductance has an alternating 
peak structure via the symmetric and anti-symmetric states because of the finite inter-dot tunnel coupling. In this 
situation, the peak heights of the conductance via the symmetric state are higher than those via the anti-symmetric 
state since the energy of the latter is further from Fermi energy than that of the former. 



3. Strong coupling regime 

Next, we consider the strong coupling regime, |eo| < \tc\- When a = 1, there is one conduction mode and the 
perfect transmission G/Gq = 1 is realized. For the parameters shown in Fig. ^1 the resonant condition (eo^ — tc^ — 
— (fi.7)^ sin^((/)/2)) is satisfied at $/$o = (2n -(- 1)/4, where n is an integer. Unlike to the weak coupling regime, there 
are no completely destructive interferences in the AB oscillations. The visibility decreases as a decreases, which is 
similar to the results for the weak coupling regime. 

It is interesting to note that the conductance is maximum at integer $/<I>o foi' weak coupling, while the conductance 
is minimum at that flux for strong coupling. This is more evident by comparing the data at a — 0.5 in Figs. I^and 
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FIG. 10: AB oscillations in the strong coupling regime, eo/h-y — 1 and tc/h-y — 
indicate a = 1, 0.5, and 0, respectively. 
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FIG. 11: AB oscillations in the double resonant levels, eo/hj = tc/hy = 0. The solid, dotted, and broken, and dash-dotted 
curves indicate a = 1, 0.9, 0.5, and 0, respectively. 

irni Such a change of the AB oscillation phase by tc bears intriguing analogy to the AB oscillation phase change by 
the energy difference between the two QDs pointed in Ref. 0. 

4.. Double resonant levels 

Here we consider a special case when eo = ic = 0, namely two energy levels of the DQD are aligned with the Fermi 
level (resonant levels). The linear conductance shows a singular flux dependence at a = 1 as discussed in Refs. |^ 
andll^l (see Fig. [TT|l . namely the linear conductance becomes zero except for $/4>o = n where G = Gg is realized. 
As in the previous sections, the visibility decreases as a decreases and finally the perfect transmission {G/Gq = 2) is 
realized at a = since two conduction modes satisfy the resonant condition (eg = 0). 

IV. TUNNELING CURRENT IN AN ASYMMETRIC AB INTERFEROMETER 



In this section, we discuss the tunneling current in an asymmetric AB interferometer, namely 7^7 ^ 7^ or (pjj ^ (pL. 
However, it is well-known that the asymmetry of the linewidth functions suppresses the Breit-Wigner resonance 
because of the asymmetric factor r'^r^/(r'^ 4- F^). Therefore, in this section, we only study the asymmetry of the 
flux. As in the previous sections, there are two conduction modes in the asymmetric systems except when a = 1 where 
only a single conduction mode contributes to the transport. The behavior of the Fano effect due to the asymmetry 
of the flux has been already investigatedS^. In their analysis of a = 1, the conductance can be expressed as the Fano 
form with flux-imbalance dependent energy scaling and Fano parameter q. We studied the effect of a less than 1 and 
found that the Fano resonance is very sensitive to the parameter a similar to the results in Sec. IIII CI (not shown). 
Here we only show the AB oscillation results with several a parameters emphasizing the difference for regimes of 
interest. 

In the situation of zero coupling {tc = 0), only the total flux influences the relative phases of the electrons tunneling 
through the two dots since there is no direct tunneling between dots. Thus, the visibility of the AB oscillations is 
independent of the asymmetry of the flux in the weak coupling limit. 
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FIG. 12: 7cr = 7l = 7, $u/^o = $/6"3?o, 'I'l/'3?o = 5$/6<l>o, = 2, and tc/h'y = —1. The solid, dotted, and broken curves 

indicate a = 1, 0.5, and 0, respectively. 
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FIG. 13: 7(7 = 7l = 7, ^u/^o = 'I'/e^o, <I>l/$o = 5$/6$o, eo/^T = 1, and tc/hj = -^^3/2. The solid, dotted, and broken 
curves indicate a = 1, 0.5, and 0, respectively. 

A. Weak coupling regime 

Here we consider the weak coupling regime, |eo| > \tc\- The AB oseillations for = 1/48 and ~ 5/48 are 
shown in Fig. ^1 When we compare Figs. l^and lT^ the period of the AB oseillations changes due to the asymmetry 
of the flux. The periods of the AB oscillations are determined by the flux threading through a smaller subcircuit. In 
Fig. ^1 the flux threading through a smaller subcircuit is $/6$0j and the period is 127r. When the ratio of the flux 
is an irrational number, the AB oscillations no longer have a well-defined period2&. For \a\ < 1, the visibility of the 
AB oscillation decreases as the \a\ decreases as well as the symmetric situation. 

B. Strong coupling regime 

Here we consider the strong coupling regime, |eo| < \tc\- The AB oscillations for = 1/48 and $l = 5/48 are 
shown in Fig. ^| When we compare Figs. llUland lT^ the period of the AB oscillations changes due to the asymmetry 
of the flux. The periods of the AB oscillations are determined by the flux threading through a smaller subcircuit 
as found in the weak coupling regime. Unlike the situation of the symmetric flux as discussed in Sec. Illl D 31 where 
completely destructive interference was absent, we could observe sharp zero conductance dips when a = 1. The 
visibility of the AB oscillation decreases as the a decreases, which is similar to the symmetric situation. 

V. CONCLUSION 

In this paper, we investigated the electron transport through a DQD system with the indirect coupling parameter 
a. We found that there arc two conduction modes except when \a\ = 1. The Fano anti-resonance is suppressed for 
|a| < 1. The visibility of AB oscillations decreases as |a| decreases and the visibility reaches zero when a — since 
the tunneling conductance is independent of the flux. Moreover, we have investigated the effects of the asymmetry 
of the flux and the tunnel coupling between reservoirs and dots. The asymmetry of the flux leads to a change in the 
period of the AB oscillations. 
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APPENDIX A: DERIVATION OF INDIRECT COUPLING PARAMETER a 

Here we show the detailed derivation of the indirect coupling parameter a discussed in Section flll Al For the model 
given in Section IlII Al the wave functions of reservoir modes and their eigenvalues are 

<I^Uk{r) = ^e^^VP sin f ^"j e-^ (Al) 



where the decay constant k is 



and the scattering phase shift satisfies sin 6* = ^2^^" ^ derived by matching wave functions C(j/) = e**^"^ + re"*'^*^ 
for y < and C,{y) = Ce"'**' for y > 0. We foimd r = — e** and C = — 2ie*^ sin|. L is the system size and d 
is the dimensionality of the reservoir, three-dimension (3D) (A:|| ~ {kx,kz) and p ~ (a;, z)) or two-dimension (2D) 
(/cy — {kx^Q) and p = (a;,0)). For the dot, we consider only the symmetric mode, namely zero angular momentum 
state, 

where cr^ + ^11^ = '^^ a-nd Dd is a constant related to the detail of the potential inside the dot and Kq is the 0th order 
modified Bessel function. The origin should be shifted to the center of each dot fi . To calculate the tunneling matrix 
element, it is convenient to obtain the Fourier transform of them, for 3D: 



and for 2D: 



Using Eq. Q for 3D system, the tunneling matrix elements are 

2/xL3/2«in(^2j ^ ^ ' 

and for 2D system are 

Therefore, using the dimensionality-dependent constant A^, we have in general 

i[;i^A,sinf^')e-^-+''=^--. (A9) 
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Using Eq. © for 3D system, the general formula of linewidtli functions at the Fermi energy /iy = ^ is 



where s — \xi — Xj\, which have a limiting form for U ^ ^jj, 

\2 A J- / r \ 3 



vUm) - [£) e-^y-y^^^'j-nik.s). (Aii) 



and for U > ^u, 



Similarly, the formula for 2D system is 

2 „fc 



and a limiting form for U ^ /if/, 



r.^(A^a) - ^ ( ) ^e-^(-+^^) J,(fc..), (AM) 



27rA2^ f L\nkf 

and for C/ > /i£/. 



We define the indirect coupling parameter as follows: 



Therefore, the parameter at the Fermi energy for 3D system with U S> 



a^{e) ^ S''^'l.,. - (A16) 



where ji is the first order spherical Bessel function and for U > /iy 



a"" {^iu) - ^ji{kFs), (A17) 



where we consider symmetric configuration: yi ~ y2 ~ Vd- For 2D system with U ^ jiui 



a^'ifiu)-- ^Mkps), (A19) 



where Ji is the first order Bessel function and for U > ^ijj 

where we consider symmetric configuration. In Ref. ITsl although similar derivation of the parameter a of 2D system 
was discussed, their result (a oc Jo{kps)) is different from ours. Our result is more general since the authors of Ref. 
ITsI ignored ky dependence of the scattering phase shift. 
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APPENDIX B: TRANSMISSION PROBABILITY 



We diagonalize the transmission matrix at the Fermi energy T(0) = G'~(0)r'^(0)G''(0)r^(0) fomid in Eq. l(TS)). 
We then obtain two eigenvalues of the transmission matrix 



Ti{e) - 
T2{e) = 



{h^f[TA{e)-TB{e)] 



(e - eo)2 - ic' - {hiY {l - «2 cos2 |] + 4(^,^)2 _ + ^os (|) 

{h^f[TA{e)+TB{e)] 



"(e - eo)2 - ic' - (^7)2 {1 - a2 cos^ (|) }] + 4(^17)2 [(e - eo) + at, cos (|) 
with the two positive functions 



2 ' 



(Bl) 
(B2) 



TA{e) = (e - eo) (1 + « cos0) + 4a(e - eo)ic cos ( | ) + (1 + a^)t/ + (?i7)'(l - a^) <{ 1 - a"" cos" 



TB{.e) = 2 



(B3) 



a{e ~ eo) cos ( - ) + tc 



Xi a- 



2(%)2(l-a2)sin2 (^0 +(e-eo)2|l-a2sin2 (^0 | + 2a(e - eo)te cos + aH,\ (B4) 



When \a\ = 1, 



rA(e)=Ts(e) = 2<i(e-eo)cos(^ , 



and the transmission probabihties are 
Ti(e) = 0, 



4(^7)2 [(e_eo)cos(|) + 



(e - eo)2 - _ (;j^)2 sin" ( |) + 4(;i7)2 (e - e^) + cos ( | 



(B5) 

(B6) 
(B7) 



Therefore, we have single conduction mode at |a| = 1. When \a\ < 1, Ta 7^ Tb and the number of the conduction 
modes is two. In fact, we have the following relation: 



Ta^ - Tb" = (1 - [l + 5(2 - + (e - _ t,)^] [l - 5(2 - 5) + (e - eo + Q^] > 0, 



(B8) 



where a = a cos 

In particular, when a = 0, the transmission probabilities are 



Tiie) = 
T2{e) = 



(e - eo - te)2 + {hj)^ ' 
(%)^ 

(e-eo + ic)2 + (;i7)'' 



(B9) 
(BIO) 



Ti and T2 represent Breit-Wigner resonances through the symmetric and anti-symmetric states, respectively. These 
are independent of the flux. 
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